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ABSTRACT 

Context. Luminous infrared galaxies (LIRGs) are an important class of objects in the low-z universe bridging the gap between normal 
spirals and the strongly interacting and starbursting ultraluminous infrared galaxies (ULIRGs). Since a large fraction of the stars in 
the Universe have been formed in these objects, LIRGs are also relevant in a high-z context. Studies of the two-dimensional physical 
properties of LIRGs are still lacking. 

Aims. We aim to understand the nature and origin of the ionization mechanisms operating in the extra-nuclear regions of LIRGs as a 
function of the interaction phase and infrared luminosity. 

Methods. This study uses optical integral field spectroscopy (IPS) data obtained with VIMOS. Our analysis is based on over 25 300 
spectra of 32 LIRGs covering all types of morphologies (isolated galaxies, interacting pairs, and advanced mergers), and the entire 
10" - 1O'"L0 infrared luminosity range. 

Results. We found strong evidence for shock ionization, with a clear trend with the dynamical status of the system. Specifically, we 
quantified the variation with interaction phase of several line ratios indicative of the excitation degree. While the [Nn]/16584/Ha ratio 
does not show any significant change, the [S n]/l/16717,6731/HQ' and [O i]/16300/H(!' ratios are higher for more advanced interaction 
stages. Velocity dispersions are higher than in normal spirals and increase with the interaction class (medians of 37, 46, and 51 km s"' 
for class 0, 1, and 2, respectively). We constrained the main mechanisms causing the ionization in the extra-nuclear regions (typically 
for distances ranging from ~0.2-2.1 kpc to ~0.9-13.2 kpc) using diagnostic diagrams. Isolated systems are mainly consistent with 
ionization caused by young stars. Large fractions of the extra-nuclear regions in interacting pairs and more advanced mergers are 
consistent with ionization caused by shocks of VsS200 km s"'. This is supported by the relation between the excitation degree and 
the velocity dispersion of the ionized gas, which we interpret as evidence for shock ionization in interacting galaxies and advanced 
mergers but not in isolated galaxies. This relation does not show any dependence with the infrared luminosity (i.e. the level of star 
formation). All this indicates that tidal forces play a key role in the origin of the ionizing shocks in the extra-nuclear regions. We 
also showed for the first time what appears to be a common log([Oi]/16300/HQ-) - log(cr) relation for the extranuclear ionized gas in 
interacting (U)LIRGs (i.e. covering the entire 10"" - 10'^ ^^ Lq luminosity range). This preliminary result needs to be investigated 
further with a larger sample of ULIRGs. 
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1. Introduction 

Luminous and ultraluminous infrared galaxies (LIRGs and 
ULIRGs) are defined as those objects with an infrared luminos- 
ity of L/R = L(8 - lOOO^/m) = 10" - lO'^Lfr, and L,r k lO'^-Lfr-, 
respectively (see lSanders & Mirabelll996tlLonsdale et al1l2006l 
for a review). They are systems which contain large amounts 
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of gas and dust (e.g. lEvans et al ] l2002h and which are undergo- 
ing an intens e star-formatio n episod e in their (circum)nuclea r 
regions (e.g. IScoville et al] |2000; Alonso-Herrero et al.ll2006l) . 
This activity is the main cause of their huge luminosity in 
about ~80% of these systems, although some contribution 
from an AGN is present and even dominant in some cases 
(e.g. Genzel et al. ' 19981: iRisaliti et all [20061: iFarrah et all 120071: 
Nardini et al 20081) . 

These systems usually present some degree of interaction 
whose importance increases with luminosity. While the majority 
of local LIRGs can be classified as isolated spirals or interact- 



2 



Ana Monreal-Ibero et al.: VLT-VIMOS integral field spectroscopy of LIRGs and ULIRGs II. 



ing pairs (e.g. Ayibas et al.l l2004t lAlonso-Herrero et al 
ISanders & Ishida"2004'). most of the ULIRGs show signs of a 
clear merging process (e^. Clements et al. 1996; Borne e t al 
20001: ICuiet all I2OOII: iBushouse et all 120021: IVeilleux et al 



tative sample of low-redshift LIRGs and ULIRG s using opti- 
cal and near-IR IFS facilities ( Arribas et al.l l2008l) . Specifically 
we used the INTE GRAL+WYFFOS facihty dArribas et at] 
119981: iBingham etTal.! 1 1994 and th e Pots dam Multi-Aperture 



2OO2I) . While (U)LIRGs are an oddity in the local Universe, Spectrograph, PMAS (lRothetal.1 l2005b in the Northern 



recent mid-infrared and submillimeter surveys show how they 
present a strong evolution with redshift, increasing their number 
by two orders of magnitude at z ~ 0.8 - 1.2 (Elbaz et al. 2002). 
Indeed they are the dominant population of the infrared selected 
galaxies at high redshift, making a signific ant contribution to the 
star-formation rate density a t 0.5 < z < 2 dPerez-Gonzalez et alj 
l20i05tlLe Floc'h et al.ll200l . 

The study of the ionization properties of the gas in these 
objects is relevant for two main reasons. On the one hand, the 
ionization is important to investigate the nature (i.e. starburst, 
AGN) of the dominant source that causes the huge luminosity 
in the infrared. In the optical, this has mainly been done via 
long-slit observ ations of the nuclear regions of large samples of 
(U)LIRGs (e.g. lKim et aD[T995Hwiieux et alJI 19991 and refer- 
ences therein). These studies established trends with the lumi- 
nosity and interaction stage, and found an increase in the fre- 
quency of AGN-dominated systems with luminosity. These re- 
sults have been recently revisited using the new optical classifi- 
catio ns provided by the use of Sloan Digital Sky Survey (SDSS) 
data dYuan et al.ll2010l) . They show that most of the (U)LIRGs 
previously classified as Low-ionization nuclear emission-line 
region (LINER), now are classified as starburst-AGN com- 
posite galaxies. The presence of an obscured AGN has been 
also revealed by the detection of ionization cones with integral 
field spectroscopy (IFS) data (e.g. Arp 299. iGarcia-Marm et alJ 
[2006). 

On the other hand, the ionization structure helps to under- 
stand how the interaction/merger process as well as the release 
of energy and material from the central source and/or starbursts 
are affecting the extended structure of the galaxies in general and 
its interstellar medium in particular. In that sense, the presence of 
Super Galactic Wi nds (SGWs) in (U)LIRGs has been suggested 
using emission ("Heckman et al.' '1990*; 'Lehnert & Heckman' 
1996[) and absorption (Heckman et al. 2000; Rupke et al. 2002, 



Hemisphere, and VIMOS dLeFevre etal.ll2003l) and SINFONI 



2005 a. b) lines. Tidally induced forces associated with the inter- 
action process itself h ave been also suggested as the cause for th e 
ionization of the gas dMcDowell et alJl2003tlColina et al.ll2005l) . 
Given the complex structure of these systems, where the selec- 
tion of a preferential direction is specially difficult, these studies 
would benefit from IFS data thanks to which it is possible to 
obtain homogeneous two- dimensional spectral information. 

Using this technique iMonreal-Ibero et aL I (l2006l hereafter 
MAC06) have studied a sample of six ULIRGs (nine galaxies), 
and found that wide areas of the extra-nuclear extended regions 
presented line rati os typical of LINERs accord ing to the diag- 
nostic diagrams of lVeilleux & Osterbrockldl987 ). In addition, it 
was shown that the velocity dispersion is positively coiTelated 
with the degree of ionization supporting the idea that shocks are 
the main cause of the ionization in these areas. However, these 
results were based on a relatively small sample, which covered 
a restricted range in luminosity (log(LjK / Lq) - 12.03 - 12.40) 
and interaction phase. 

In this paper we extend that study to a larger sample of 32 
systems, which cover the entire \og(LiRl Lq) - 1 1.00 - 12.00 lu- 
minosity range (i.e. the LIRGs range), and the different interac- 
tion types (i.e., isolated galaxies, interaction pairs, and mergers 
remnants). 

The present study is part of a wider project devoted to the 
study of the internal structure and kinematics of a represen- 



(Eisenha uer et al.l l2003') in the southern one. The coiTesponding 
catalogs for t he PMAS, INTEGRAL and V IMOS samples can 
be fou nd in lAlonso-Herrero et aL i d2009l) . |G arcia-Marfn et al.l 
d2009 ') and Rodriguez -Zaurm et al. (in prep.), respectively. 

The paper is structured as follows: in Sect. |2]we describe 
the sample used in this work as well as the characteristics of 
the instrumental configuration and technical details regarding 
data reduction and analysis; Sect.|3]quantifies how the ionization 
degree varies with interaction stage and constrains the possible 
mechanisms that cause the ionization of the gas. Finally, a com- 
parison with the previous results for ULIRGs and a discussion 
about the origin of the ionization produced by shocks in terms 
of the star formation and the interaction process are presented. 

Throughout the paper, a cosmology with 70 km s"' Mpc"', 
Q.M - 0.3 and Qa - 0.7 is assumed. 



2. The data 

2.1. The sample 

The present sample is drawn from the VI MOS IFS sample of 
(U)LIRGs presented in lArribas et alj d2008 , hereafter. Paper I). 
Specifically, it includes all the LIRGs listed in that paper, except 
for IRAS F10173-I-0828, for which no emission lines were de- 
tected. That implies 32 systems, of which 13 were isolated, 11 
interacting systems (9 pairs and 2 triple), and 8 advanced merg- 
ers (i.e. classes 0, 1, and 2, respectively, according to Paper I 
terminology). The methodology followed to perform the mor- 
phological classification is described below. 

Figure[T]summarizes the distribution in luminosity, distance, 
and interaction type of the systems in our final sample. The mean 
distance of the sample is 87 Mpc, which leads to a mean lin- 
ear scale of ~400 pc arcsec"'. Taking into account the VIMOS 
spaxel size of 0'.'67, this translates into a mean linear spatial sam- 
pling of the source of ~270 pc spaxel We refer the reader to 
Table 1 in Paper I for the basic properties of the individual sys- 
tems. 



2.2. Morphological classification 

Our morphological/ merging classification is a simplified version 
of that proposed bv lVeilleux et al.l d2002h for ULIRGs, who di- 
vided their sample in five classes (and four sub-classes). We only 
considered three main classes (i.e. "isolated": 0, "interacting": 1, 
"merger": 2) to reduce uncertainties associated with sorting. In 
particular, the different morphological classes considered are 

- Class 0: Objects that appear to be single isolated galaxies, 
with a relatively symmetric morphology and without evi- 
dence for strong past or ongoing interaction. 

- Class 1 : Objects in a pre-coalescence phase with two well 
differentiated nuclei separated by a projected distance of 
D>1.5 kpc. For those objects classified as 1, it is still possible 
to identify the individual merging galaxies and their corre- 
sponding tidal structures due to the interaction. The limit of 
1.5 kpc was chosen because theoretical models predict a fast 
coalescence phase after the nuclei become closer than t hat 
distance (e.g. lMihos & Hernquistl[T996tlNaab et all2006h . 
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Fig. 1. Histograms showing the luminosity (left) and distance (right) distributions of the systems in the sample. Light, intermediate 
and dark gray areas indicate the total number of isolated galaxies (class 0), interacting systems (class 1), and merger remnants (class 
2), respectively (see text for a more detailed definition of the different morphological types). 



- Class 2: Objects with a relatively asymmetric morphology 
suggesting a post-coalescence merging phase. They may 
have two nuclei separated a projected distance of D< 1 .5 kpc. 
For objects classified as 2, it is not possible to identify the in- 
teracting galaxies individually. 



For simplicity we will refer to objects of classes 0, 1, and 2 
as "isolated", "interacting", and "merger" systems. 

The classification was primarily based on the Digital Sky 
Survey (DSS) images, which are available for all the sources 
on the NASA Extragalactic database (NED). In most of the 
cases these images are sufficient to classify the objects of the 
present sample. However, for seven objects (less than 25% of 
the sample), the DSS images were not conclusive. In these cases 
(IRAS 08355-4944, IRAS F08520-6850, IRAS F10038-3338, 
IRAS F12116-5610, IRAS F13001-2339, IRAS F17138-1017 
and IRAS F21453-3511) high-resolution HST images were 
used as supplementary information, mainly from program 
ID: 10582 (IP: Evans). For another six objects the HST images 
just confirm the morphological classification derived using the 
DSS images. In Figure [T4l (on-line) we present the images used 
for the classification. In some cases the dynamic range is large 
and the different relevant features appear at quite different inten- 
sity levels, so we recommend the reader interested in a particular 
case to directly download and display the images. 

In order to estimate the uncertainty associated with this 
classification process, we thrice classified the whole sample 
independently. The level of agreement was higher than 90%. 
However, there were intrinsically difficult cases for which we 
agreed they have an uncertain classification. These cases are: 
IRAS 09437-H0317 (IC 563/IC 564) (1/0), IRAS F121 16-5615 
(2/0), IRAS F13001-2339 (2/1), and IRAS F17138-1017 (2/0), 
and they will not be considered later when analyzing the global 
behavior of the different classes. 

Further details on the morphology of these objects may be 
found in Rodriguez-Zaurm et al. (in prep.). 



2.3. Observations and instrumental set-up 

A detailed description of the observations is provided in Paper I. 
Here we briefly recall the observational set-up and characteris- 
tics of the data. 

The data were obtained in service mode with the VIMOS- 
IFU at VLT during semesters 76, 78, and 8 1 . We used the HR- 
Orange configuration, which covers the 5250 - 7400 A spectral 
range with a resolution of 3 400. It provides a field-of-view of 
27" X 27" with 0'.'67 per spatial element (spaxel) making a total 
of 1 600 spectra per pointing. Each galaxy was observed using a 
4-pointing dither pattern with a relative off'-set of 2"7 (i.e. four 
spaxels) to minimize the eff'ect of dead fibers and thus, providing 
an effective field-of-view of about 29'.'7 x 29'.'7. Details about 
data reduction and line fitting can be found in Paper I. 

2.4. Emission iine data 

We used in our analysis the [Oi]^6300, Ha, [Nii]ii6548,6584, 
and [S ii]/l/l6717,6730 emission lines, which were fitted to 
a single Gaussian component. For some galaxies, small ar- 
eas of the Ha emission showed evidence for two or more 
kinematically distinct components. These particular areas of 
double components are generally associated with nuclear re- 
gions (i.e. IRAS F08520-6850, IRAS Fl 3229-2934, the east- 
ern member of IRAS F14544-4255 and IRAS 21453-3511) 
which are not used in the present analysis (see below). 
In some cases these also affect small areas in the extra- 
nuclear regions (IRAS F06592-6313, IRAS F07160-6215, 
IRAS F10409-4556, IRAS F13229-2934, IRAS 08424-3130, 
IC 564, the western member of IRAS F14544-4255, 
the northern and central members of IRAS 18093-5744, 
IRAS F04315-0840, IRAS 10257-4338, IRAS 17138-1017 
and IRAS 21453-3511). For a given galaxy, these regions rep- 
resent typically less than 5% of the data. In those cases where 
the two components were clearly distinguishable, the dominant 
component was used in the present analysis, while in those of 
strong blending we used the results from the one-profile fit. A 
more detailed analysis of the regions where multiple components 
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are identified is out of the scope of the present paper and will be 
presented with the full kinematic analysis in a future paper 

Thus a measure of the flux in the diff'erent emission lines as 
well as one independent measurement of the central wavelength 
and FWHM (full width at half maximum) is obtained for each 
spectrum. The FWHM has been corrected from the instrumen- 
tal width measured with the strong [Oi]/i6300 A sky line, and 
then translated into velocity dispersions. For the VIMOS con- 
figuration used, the spectral resolution translates into an instru- 
mental width cr,„j ~ 38 + 8 km s ^ The corresponding values 
for the weaker [O i]/l6363 A sky line are cr,,,, ~ 33 + 10 km s"' . 
Therefore, the profiles of the emission lines will be considered as 
resolved if their observed profile is wider than 50 km s"' (ctq/,,). 
In general the [Oi]/I6300, [S ii]^^6717,6731, [Nii]^^6584, and 
Ha emission lines define a sequence of increasing S/N and, 
therefore, the [N ii]/16584/Hq' maps cover a larger area than 
those generated with the other line ratios. While uncertainties 
in the line flux for the strongest emission lines (Ho-) in high 
surface-brightness regions are about 10%, weaker lines like 
[S ii]/l/l67 17,6731 and [Oi]/I6300 can have large uncertainties 
due to the lower S/N. However, no lines with flux uncertainties 
larger than 30% are used in the analysis presented in this paper. 
So, the typical uncertainties in the line ratios varies from about 
15% for high S/N lines to about 30% for low surface brightness 
and when a weak line is involved. 

Similarly to MAC06, we only used the extra-nuclear regions 
of our galaxies for the analysis in this paper. To do so, we ex- 
cluded the nuclear regions, defined as those within the central 
3 spaxelx3 spaxel (2'.'0 x 2'.'0) where, in addition to the complex 
line profiles mentioned above, contamination due to a dust en- 
shrouded AGN could aff'ect the results. For the typical distances 
of our sample of LIRGs, this corresponds to a region of about 
~800 pcx800 pc in size. For the particular case of IRAS F043 15- 
0840, where the two galaxies are coalescing, the area with emis- 
sion lines presenting double components was a bit larger than the 
standard 3 spaxelx3 spaxel one. Because here the one Gaussian 
approach did not properly reproduce the line profiles, we we 
made an exception and also masked the corresponding spaxels 
(see Fig.|2]i. 

The final total number of available data points, mean and 
standard deviation for each individual pointing for the three line 
ratio used as well as the velocity dispersion inferred from Ha are 
shown in Table [1] 

3. Results and discussion 

3.1. General characteristics of tiie excitation maps 

Figure E] presents the [Nu]A6584-/Ha, [S ii]^^6717,6731/Ha', 
and [Oi]/16300/Hq' line ratio maps for each individual 
VIMOS pointing. Although these maps onl y represent one 
axis of the classical BP T diagrams (Bald win et al.l ll98U 
IVeilleux & Osterbrockl ll987). they are good tracers of the exci- 
tation when studying the general properties of the extra-nuclear 
ionized regions of LIRGs. In fact, in flie X-[Oiii].i5007/H/? di- 
agnostic diagrams (where X is any of the above line ratios), the 
boundaries distinguishing regions mainly photo-ionized by stars 
(i.e. H ii-like) and those of LINER type are nearly vertical for the 
expected values of log([Oiii]/l5007/Hy6) in these extra-nuclear 
regions (~ -0.2, see below). Figure 2 also shows the Ha sur- 
face brightness maps. These maps are presented here mainly 
as a reference for the excitation maps, and they will be dis- 
cussed in detail elsewhere (Rodriguez-Zaurm et al. in prep.). The 
different area covered by the excitation maps is a direct con- 



sequence of the diff'erent signal-to-noise (S/N) of the involved 
lines. The [Nii]/16584/Hq' map encompasses the largest area 
in all the cases, reaching areas of low Ha surface brightness 
where no good S/N data are available for the [Oi]/I6300 and 
[Sii]/l/l6717,6731 lines. 

The star-forming regions of relatively low excitation are 
well identified in the maps. These are generally high Ha 
surface brightness regions associated with large scale struc- 
tures like spiral arms (e.g. IRAS Fl 3229-2934), rings (e.g. 
IRAS Fl 1255-4120), or tidal tails (e.g. IRAS 06076-2139). 
However, knots of star formation are also found in isolated exter- 
nal regions (e.g. IRAS F121 15-4656, IRAS F07027-6011 S). 
In some cases these external regions define a chain suggestive 
of star formation along tidal tails (e.g. IRAS Fl 1506-3851, 
IRAS 06076-2139, Paper I). 

Galaxies in double or triple systems may have quite differ- 
ent excitation conditions, indicating that the interstellar medium 
of the individual objects have different properties. For instance, 
while the external regions of the central and southern component 
of IRAS F06259-4708 are dominated by star formation, for the 
northern galaxy these regions show higher excitation (similarly 
for IRAS F06076-2 139). 

Interestingly the external regions associated to diffuse low 
surface brightness Ha emission have relatively high excitation. 
In these regions the line ratios could, in principle, be affected 
by an underlying absorption spectrum. I ndeed, preUminary re- 
sults for our PMAS sample of LIRGs (lAlonso-Herrero et al.l 

2009) show that for spectra with EW(Ha)<lO A, the fine ra- 
tio can decrease ~0.2 dex once the contrib ution of an old stellar 
population has been taken into account (iAlonso-Herrero et al.l 

2010) . In order to quantify the effect of a possible compo- 
nent in absorption, we created maps including a correction 
of EWabs = 2 A. [Sii]^^67 17,673 l/Ha and [Oi]/16300/Hq' 
line ratios are affected by about 0.2 dex for only a handful 
(<2%) of spaxels, and by less than 0.1 dex, for most of the 
data. For the [N ii]/16584/Hq' line ratio, this effect could reach 
~ 0.3 dex in the very low surface brightness areas, but does 
not change in any case the observed ionization structure that 
we describe here. These regions of low Ha surface brightness 
and relatively high excitation are found in all kind of objects: 
structured systems with rings (IRAS 10567-4310) and spiral 
(IRAS F04315-0840) morphologies, highly disturbed systems 
(the central member of IRAS F07 160-5744), and galaxies with 
relatively round/ elliptical shape (IRAS Fl 1506-3851). In ob- 
jects with round morphologies this transition from low to high 
excitation translates into the appearance of external rings of 
high excitation (IRAS Fl 1506-3851, IRAS F121 15-4656). The 
transition from low to high excitation is also observed in the 
triple system IRAS 18093-5744 where an extension (tidal tail?) 
from the main body of the central galaxy towards the northern 
galaxy is identified. In some cases the high-excitation regions are 
found close to the nucleus and/or along preferential directions, 
and defining structures suggestive of cones. Some examples are 
IRAS F10409-4556, ESQ 297-G012 or the southern member 
of IRAS 12042-3140. 

For the rest of the paper we will focus on the study of the 
extra-nuclear excitation conditions using the presented Une ra- 
tios as well as the velocity dispersions. 
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Fig. 2. Maps for Ha and [Nn]^6584/HQ' [S ii]^^6717,6731/Ha', and [Oil^eSOO/Ha line ratios for our LIRGs sample. Each row 
shows a VIMOS pointing. The nuclear masked areas not considered in this work are represented in the line ratio maps as white 
squares. Ha intensity is displayed in arbitrary units and using a logarithmic scale to better emphasize all the morphological features. 
The physical scale of 2 kpc at the distance of the galaxy is represented by the straight line in the bottom right corner Axes scales 
are in arcsec and orientation is as usual: north up, east to the left. 
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3.2. Dependence of the excitation conditions of tlie 
extra-nuclear ionized gas witli tlie interaction class 

In order to investigate possible changes in the excitation prop- 
erties of the ionized gas as a consequence of interactions and 
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Fig.|2l continued. 



mergers, the distribution of the emission line ratios for the entire The [Nn]/I6584/Ha, [S n]/l/l6717,6731/Ha, and [Oi]/16300/Hq' 
sample of LIRGs was divided into three main groups according emission line ratio distributions for these classes of galaxies as 
to their morphology; isolated, interacting, and merger remnants, well as for the entke sample are shown in Fig. [3] Galaxies with 
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Fig.|2l continued. 



dubious classification (see Table [U were excluded from the his- tions are indicated in Table |2] We also indicate in each panel 
tograms or the individual groups. The number of data points the boundaries betw een ionization c aused by stars and by other 
as well as the median and standard deviation of these distribu- mechanisms in the iBaldwin et al.l (ll98lL hereafter BPT) dia- 
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grams according to the b oundaries proposed by differe nt authors. H ii regions, and nuclear (or integrated) data for starburst galax- 
The original boundaries dVeilleux & OsterbrocU 1987 1) were em- ies and different kinds of active galaxies (i.e. Seyfert 2, LINERs, 
pirically determined using a sample that includes extragalactic narrow-line radio galaxies and what would be called today "in- 
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termediate objects"). Despite a recent new discussion about the ical boundaries for the three BPT diagnostic diagrams ii) they 
location of these boundaries, we considered that they provide facilitate a possible comparison of the line ratios presented here 
valuable information because i) they are the only set of empir- with those in previous works. The other complete set of three 
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boundaries is the one proposed bv lKewlev et al.l (12001b'). Using 
a combination of photo-ionization and stellar populations syn- 
thesis models they determined the extreme cases under which 
line ratios can be explained via photo-ionization caused by stars. 
Recently, the empirical borders associated to the [Nii]/16584/Ha 
line ratio have been up-dated by using the Sloan Digi tal Sky 
Survey (hereafter SDSS,.Ys rkeLal. 2000) data ( KaufFm ann et al.l 
120031: IStasinska et al.ll200^ There is no similar up-date for the 
diagnostic diagrams involving the [S n]/l/1671 7,673 l/Ho' and the 
[Oi]/16300/Hq' line ratios. 

In order to locate the ionization type boundaries in Fig, [3] 
we assumed a log([Oiii]/l5007/Hy6)=-0.2. This is the median 
observed value for the samp le of LIRGs sample presented in 
lAlonso-Herrero et al.l (l200i ). Individual values range between 
-0.6 and O.lQ (excluding NGC 7469, which is known to host 
a luminous Seyfert 1 nucleus that largely affects its integrated 
spectrum). As is illustrated in Fig. |3] with gray bands, the bound- 
aries change very little in this [Oiii]/t500 7/H8 range (typically 
~0.1 dex for the empirical boundaries of Veilleux & Osterbrock 
([1987 ). Stasiriska et al. (2006) and Kauffmann et_.al. ( 2003) and 
~0.2 dex for the theoretical boundaries of Kewl ev et al.l (1200 fbb ) 
making this log([0 iii]/l5007/Hj6)=-0.2, a reasonable assumption 
for establishing the mean borders. 

The line ratio distributions for the whole sample (Fig. [3] 
upper row) are not symmetrical around a mean value, but 



' These ratios have not been corrected for the underlying stellar 
absorption in Hp. Such an absorption could t ypically decrease the 
[Om]i5007/H/3 line ratio by -0.0-0.2 dex (see lAlonso-Herrero et al] 
120091 for details). 



show a significant wing towards high values. This suggests 
that a fraction of regions have a relatively high ionization (i.e. 
LINER-like). According to the Veilleux & Osterbrock bound- 
aries, the percentage of data presenting a LINER spectrum cor- 
responds to 19%, 3 1 %, and 35% when using the [N u]A65M/Ha, 
[Sn]/l/167 17,673 1/Hq', and [Oi]/i6300/Ha, respectively. Using 
the boundaries proposed by Kewley et al.^ ( .2001b) , the num- 
ber of spaxels that cannot be explained as purely ionized 
by stars is significantly smaller (4%, 2%, and 25%, for the 
[Nii]^6584/Hq', [Sii]/l/l67 17,673 1/Ha, and [Oi]^6300/Ha line 
ratios, respectively). The significant differences in the percent- 
ages when comparing the three line ratios could be due to 
the fact that these distributions do not come from exactly the 
same regions (i.e. set of spectra). Indeed, the [Oi]/I6300/Ha' and 
[S ii]/1/16717,6731/Hq' data points are restricted to a smaller re- 
gion than the [Nii]/16584/Hq' ones because the [Oi] and [Sii] 
lines have on average lower S/N than the [Nii] line. However, 
when these distributions are generated with data points from 
the same regions for the three lines (not shown), we find 
similar differences between the percentages when using the 
Kewley et al. boundaries (1%, 1% and 13%) and even larger 
ones when using the Veilleux & Osterbrock boundaries (7%, 
21% and 22% for the [Nii]/16584/Ha', [S ii]/l/l67 17,673 1/Ha, 
and [Oi]/16300/Hq' line ratios, respectively). This confirms that 
the [Oi]/16300/Hq' distribution has a higher percentage of data 
points with high ionization than the other two distributions. That 
the [S ii]/l/l6717,6731 and [Oi]/I6300 emission is enhanced by 
shocks and therefore the [S ii]/l/l67 17,673 l/Ha and, especially, 
[Oi]/I6300/Ho' ratio are better tracers of shock-induced ioniza- 
tion (e.g. Dopita & Sutherland 1995) suggests a significant pres- 
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Fig. 3. Distributions of the different LIRG subsamples according [Nn]A65S4/Ha, [S ii].1.16717,6731/Ha, and [Oi]/16300/Hq' (first, 
second and third columns, respectively). Distributions are presented for the entire LIRG sample (1" row), for class / isolated 
galaxies (2"'' row), class 1 / interacting (3"' row), and class 2 / advanced mergers (4''' row). The different vertical lines indicate the 
boundaries between ionization by stars and b y other mechanisms proposed i n different wor ks and assuming a log([Om]/l5007/Hj6) 
- -0.2. The code used is: continuous line - Veilleux & Osterbrock ('1987); dotted line - iKauffmann et alJ 12003); dashed line - 
iKewley et al. ( 2001b); three dots/dashed line - Stasihska et al. (2006). The gray bands show the variation of a given boundary when 
the assumed [Om]/I5007/Hy6 value changes the full range observed in the Alonso-HeiTero et al. LIRGs sample (i.e. -0.6,0.1). 



ence of this type of ionization in the extra-nuclear extended re- 
gions of LIRGs. This will be explored in more detail below. 

A clear distinction emerges in the excitation properties, i.e. 
in the distributions of the emission line ratios, when LIRGs are 
separated according to their interaction class (Fig. [3]). While 
the [Nii]/16584/Hq' line ratio distributions are similar with a 
median ratio of ~ -0.37 + 0.18 dex, typical of Hii, the 
[S ii]/l/l67 17,673 1/Ha line ratio shows on average higher values, 
i.e. higher excitation, in the more dynamically perturbed sys- 
tems (class 1 and 2). Indeed, median(±standard deviation) val- 
ues of -0.55 + 0.13, -0.49 + 0.16, and -0.48 + 0. 17 are found for 
classes 0, 1, and 2, respectively. Similar values are found after 



assuming a coiTection for a Ha component in absorption with an 
equivalent width of 2A. 

A more extreme change is evident in the [Oi]/I6300/Ha' 
ratio (i.e. shock tracer) with median values of -1.56 + 0.28, 
-1.43 ± 0.33, and -1.35 ± 0.35 for classes 0, 1, and 2, respec- 
tively. In addition, LIRGs belonging to both classes 1 and 2 
seem to have a double-peak distribution, with values around - 
1.4 (i.e. Hii-like) and -0.9 (LINER-like). In all the cases, the 
Kolmogorov-Smirnoff test allows us to reject the possibility that 
these distributions come from the same parent distribution even 
in the case of the [N ii]/16584/Hq' line ratio. 
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Fig. 4. Percentage of data for th e different interaction classes 
classified as LINERs using the Veilleux & Osterbrock' (11987b 
line ratio boundaries (up) and Kewlev et al. (2 001bl) ones 
(down) and assuming [Oiii]/I5007/Hy6=-0.2. Measured values 
are indicated with crosses, diamonds and triangles for the 
[Nii]/16584/Hq', [S ii]^^67 17,673 1/Ha, and [Oi]^6300/Hq' line 
ratios respectively. Data corresponding to the different interac- 
tion class have been joined with dashed, dotted and continuous 
lines, respectively for an easier reading of the figure. 
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Fig. 5. Variation of the median of the average line ratio per 
galaxy with the interaction class. We use the same line and color 
code as in Fig.|4]to distinguish among the emission line ratios. 



Moreover, when using the Veilleux & Osterbrock boundaries 
the percentage of data points with [S ii]/l/l6717,6731/Ha' and 
[Oi]/16300/Hq' line ratios in the LINER range increases by fac- 
tors of two and three with the interaction class with respect to the 



isolated galaxies (class 0), representing about -17% and -33% 
of the regions in class 1 and 2 galaxies, respectively. When us- 
ing the Kewley et al. boundaries, only the percentage for the 
[Oi]/l6300/HQf line ratio shows a significant increase (see Fig. 
0. 

As shown in Table [T] the number of data points associated 
with a galaxy can range from a few tens to a few thousands, 
depending on the line ratio. Given this range of two orders of 
magnitude, one might wonder if the distributions presented in 
Fig. |3]are biased due to the contribution of a few galaxies. To 
check this possibility, we generated the distributions for the three 
line ratios and the different morphological groups, but without 
the pointing with the largest number of spaxels. These were 
IRAS F22132-3705 for class 0, ESO 297-GOl 1 for class 1, and 
IRAS 10257-4338 for class 2. In this way, the remaining point- 
ings contribute in a similar manner to the distribution (i.e. with 
several hundreds of data points for the [Nii]/l6584/HQf line ra- 
tio in most of the cases and with many tens - a few hundreds 
for the [Sn]/l/l671 7,673 1/Ha and [Oi]^6300/Hq' line ratios). 
The distributions (not shown here) were similar to those for the 
whole sample, with differences between the medians ~0.05 dex. 
Moreover, we checked how the median of the average line ratios 
per galaxy changes with the interaction class. This is shown in 
Fig.|5] In this way, every pointing contributes with one data point 
per line ratio. As was found for the data distributions, this figure 
shows how, while the [N ii]/16584/Hq' line ratio remains constant 
with the interaction class, the [S n]/l/167 17,673 l/Hor and, espe- 
cially, the [Oi]/i6300/HQ' line ratio increases with the degree of 
interaction. 

Although it would have been interesting to compare the 
present line ratio distributions with those obtained for other 
kinds of galaxies (e.g. normal spirals, less powerful starbursts), 
there is a clear lack of data comparable to those presented here. 
In particular, a meaningful comparison requires galaxies with 
line ratios measured at similar spatial scales as those sampled 
here. For example, line ratios derived for the SDSS0, the largest 
extragalactic survey up to date, were derived from spectra ob- 
tained through 3-arcsec diameter fibers. For galaxies at similar 
redshift like the present sample, this would be representative of 
the ionization in the nuclear regions, and thus, not directly com- 
parable with our analysis on the external areas. For galaxies at 
larger distances - say, 0.04 < z < 0.10 - the spatial sampling in 
the SLOAN galaxies would be much larger (Q) ~ 2.4 - 5.4 kpc) 
and thus again not directly comparable. The only galaxy with 
published data comparable to those presente d here is the spira l 
NGC 5194 (also known as Whirlpool Ga/axv. iBlanc et al.l2009l) . 
a spiral galaxy interacting with a dwarf galaxy. At its distance 
(~8 Mpc) and with the size of the VIRUS-P fibers (© ~ 4'.'3) the 
spatial sampling is comparable with that for our sample within a 
factor of ~3. Figure|6]shows the distribution for the two available 
line ratios for this galaxy. While the [N ii]/16584/Hq' distributions 
for the LIRGs and NGC 5194 are similar, LIRGs have larger 
[S ii]/1/16717,6731/Hq' ratios. Because this comparison was done 
only with one galaxy, this result needs to be revised when line 
ratios in the external areas of larger samples of spiral galaxies at 
similar linear spatial resolution become available. 

3.3. Anti-correlation between excitation and Ha luminosity 

Maps presented in Sect. 13.11 (Fig. |2]i display an increase on 
the line ratio in those areas with low Ha surface brightness. 
Similar results have been observed in other environments, like 



See http : //www . mpa-garching . mpg . de/SDSS/. 
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Fig. 6. Same as Fig. |3]for NGC 5194 (only [Nii]^6584/Hq' and [Sii]/1^6717,6731/Hq' line ratios). 




Fig. 7. Relation between the Ha surface brightness and the [Nn]A6584/Har (left), [S n]/l/167 17,673 1/Ha (middle), and 
[Oi]/I6300/Ha' (right). Units are in an arbitrary scaling. The dash ed (continuous) verti cal bars represent the covered line ratio 
ranges for the sample of galactic H ii regions (DIG areas) studied bv Madse n et al.l (12006'). Note that because its [S ii]/Hq' line ratio 
only includes the A61\l sulfur line, we added an offset of 0.25 dex to their values to take into account the contribution from the 
/1673 1 line. Also the line ratio range s for the eDlG areas in a sample of nine edge-on spiral galaxies is shown with dotted hne for 
comparison dMiller & Veilleuxll2003h . Dotted-dashed lines represent our estimated observational limits. 



our Galaxy (iMadsen et al.ll2006 



Revriolds et al.lll999h. or other 



spiral galaxies (IColhns & Randll200ll: iMiller & Veilleuxll200l 
iBlancet all 12001 " in the so-called diffuse ionized gas (DIG) 
or warm interstellar medium (WIM). In this section we fur- 
ther explore this result by looking at the relation between the 
ionization degree and surface brightness for the different in- 
teraction groups. This is shown in Fig. |7] where each data 
point represents the information from an individual spaxel. 



As can be seen in the figure, while the [Nn]/(6584/HQ' cov- 
ers almost three orders of magnitude in Ha surface bright- 
ness, the [Sii] ^^671 7,673 1/Ha and the [Oi]/16300/Hq' line ra- 
tios are restricted to only about two. Also, low values of 
[S n]/l/167 17,673 1/Ha and [Oi]/16300/Ha are only found at high 
Ha surface brightness. We estimated our observational limit by 
looking at the typical values and uncertainties measured for our 
[Nii]/16584/Ha line ratio. Then we allowed for a maximum un- 
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certainty in the logarithm of the line ratios of 0.4 and assumed 
that the S/N scales with the root square of the signal. The de- 
rived observational limit is shown in Fig. Q which shows how 
these effects are caused by the sensitivity limit of the data, as 
was pointed out in Sect. 13.11 That is, the relatively low S/N of 
the [S ii]/i/l6717,6731 and [Oi]/I6300 emission lines in the outer 
parts prevent us from measuring low [S ii]/1/16717,6731/Hq' and 
[Oi]/16300/Hq' line ratios at low Ha surface brightness. Figure 
|7] shows that independently of the line ratio, high values are 
found at low surface brightness as it happened in our Galaxy 
and other spiral galaxies. Moreover, typical distances where the 
high line ratios are found (from ~400 pc up to ~6 kpc) are com- 
parable with the distances where this anti-correlation between 
excitation and surface brightness has been fou nd (up to ~8 kpc 
ICoUins & Randll200lllMiller & Veilleuxll2003l) . 

The figure also shows the ranges of [Nn]/16584/Ha and 
[S ii]/l/l6717,6731/Ha' line ratios measure d in a sample of H n 
regions and DIGs areas in our Galaxy by Madsen et al.l (|2006h . 
as well as for extragalactic DIG (Miller & Veilleux 2003|). Part 
of the [Nii]/l6584/HQf and [S ii]^^67 17,673 l/Ha line ratios are 
compatible with those expected for Hn regions. However, there 
is a large number of spaxels with line ratios similar to what it is 
observed in the DIGs areas. 

H II regions are generally understood as photo-ionized by the 
young stellar populations within them. On the other hand, DIGs 
seem more d ifficult to explain only via photo-ionization (e.g. 
iMiller & Veil leux 2003). If they are interpreted only as photo- 
ioni zation areas, there is a need for extra heating fRe vnolds et al.l 
11999; Mathis 2000) and even with this extra heating, it is dif- 
ficult to reach [N ii]/16584/Hq'>1 fine ratios. Another possibil- 
ity would be that these line ratios are the composite effect of 
photo-ionization and shocks or turbulent mixi ng layers (TML, 
IColUns & Ra"n d 2001 ; ' Miller & Veille ux"2003'). In the next sec- 
tion we will explore the role of these different ionization mech- 
anisms by comparing our measured line ratios with the predic- 
tions of the models. 



3.4. Ionization meclianisms in fhe extra-nuclear ionized 
regions: young stars and shocks 

To investigate the nature of the ionization mechanisms present 
in the extended extra-nuclear regions, the [S ii]/l/l6717,6731/Ha' 
vs. [Nii]^6584/Ha'and [Oi]^6300/Hq' vs. [N ii] ^65 84/Ha' diag- 
nostic diagrams for (thousands) of data points sampling these 
regions were compared with model predictions for H n regions, 
shocks and AGNs. State-o f-the-art models for evolvi ng Hn 
regions (Dopita et alJ l2006l) . shocks ( Allen et al.' '2008*), TML 
(Slavin et al. 1993) and dusty AGNs (Groves et al. 2004) cov- 
ering a wide range of physical parameters (ionization, density, 
stellar ages, shock velocity, etc.) and metallicities were consid- 
ered. The results for the galaxies were separated according to 
their interaction class, and set of models best representing the 
data are given in Figs.[8]and|9]for the [S ii]/l/167 17,673 l/Ho' vs. 
[Nii]^6584/HQ'and [Oi]^6300/Hq' vs. [Nii]^6584/Hq' diagnos- 
tic diagrams, respectively. Galaxies with dubious classification 
(see Table[T]l were excluded. The data associated with each indi- 
vidual pointing have been attached as On-line Material. 

From these diagrams it is clear that although the ionization 
properties of most of the extended, extra-nuclear regions are con- 
sistent with those of H ii regions, there is a trend with the interac- 
tion class. As the interaction progresses (e.g. moving from class 
to class 2), a larger fraction of extra-nuclear regions shows a 
clear shift towards an increased excitation, as traced by the in- 



crement of the [Sii]A/16717,6731 and, in particular, [Oi]/16300 
emission with respect to Ha (see Figs.[8]and|9]l. 

Although individual galaxies within an interaction class or 
regions within galaxies can behave differently, a direct compari- 
son of all the available data (represented by thousands of points 
in the diagrams of Figs.[8]and|9]i with the model predictions al- 
low us to obtain some general conclusions about the ionization 
mechanisms playing the main role in the extra-nuclear regions 
of LIRGs. 

The first conclusion that can be drawn from a compar- 
ison of Figs. [8] and |9] is that TML are not playing a ma- 
jor role in the ionization of the extra-nuclear regions of 
this sample. Models can marginally reproduce the observed 
[Nii]^6584/Hq' vs. [S ii]^^6717,6731/Ha', but clearly under- 
predict the [O i]A6300/Ha line ratio for a given [N n]A6584/Ha. 

Instead, most of the regions are consistent with ionization 
by young stellar populations with ages of 3 Myr or less and 
metallicities twice solar. The measured Ha emission line equiv- 
alent widths in these regions are in the 40 to 300 A range 
(Rodrnguez-Zaurm et al. in prep.), and therefore do not cor- 
respond with the strong Ha lines (equivalent widths of 1000 
to 2000 A) expecte d in very young (< 3 Myr) stellar clusters 
(STARBURST99, Leitherer et aljfT999 ). This apparent discrep- 
ancy is like ly due to contamination by evolved stellar popu- 
lations (see lAlonso-Herrero et alJ[2010l) . Because the physical 
scale of each resolution element (determined by a seeing of 
~ r.'2) would be about 480 pc, the measured equivalent widths 
are likely the result of different stellar populations where the 
old stars present in these regions contribute to the optical con- 
tinuum, but not to the emission lines. Moreov er, typical sizes 
for giant Hii regions would be <400 pc (e.g. Kennicuttl ll984l : 
'Alonso-Herrero et al. 112002). Because our typical resolution is 
slightly higher, a single spaxel can sample both pure DIG and 
regions of very young star formation. This implies that mea- 
sured line ratios, specially the [Nii]/16584/Ha one, are higher 
than those for pure H ii regions and thus a comparison with mod- 
els suggests too young stellar populations. 

For regions showing a log([S ii]/l/i67 17,673 1/Ha)> - 0.6, 
the model predictions considering ionization mechanisms dif- 
ferent from the TML's ones overlap and therefore no firm con- 
clusions can be made based on the [S ii]/l/l6717,6731/Ha vs. 
[Nii]/16584/Ha alone (see Fig. [8]). However, much of this over- 
lap disappears when the [Oi]/I6300/Ha ratio is considered. In 
particular regions showing log([0 i]/l6300/Ha)> -1.6 are best 
explained by the presence of shocks with velocities of less than 
200 km s"^ and metallicities 1-2 solar (see Fig.|9]l. 

We note that a detailed comparison of our data with the 
models in Figs. [8] and |9] may lead to some inconsistencies. 
For instance, the [S ii]^^6717,6731/Ha- [N ii]/l6584/Ha dia- 
gram presents less data points in the area of shocks of Vj ~ 
200 - 300 km s"' than the [Oi]^6300/Ha- [Nii]^6584/Ha 
diagram. These disagreements may be due to the intrinsic 
difficulty in modeling so me emission (e.g. [Oi]/i6300 line 
Dopita & Sutherland! 1995h as well as possible observational ef- 
fects. The areas sampled by a given spaxel include both star- 
forming regions and diffuse ionized gas, which affect the line 
ratios differently. In any case, these figures should be read look- 
ing for general behaviors rather than for detailed comparisons. 

Taking this into account, we considered the [Oi]/i6300/Ha- 
[Nii]/16584/Ha diagram as our main indicator to establish the 
most important ionization mechanism for individual galaxies in 
Table [3] In some cases this was complemented with other indi- 
cators, especially for the galaxies with no [O i]/l6300 detections. 
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Fig. 8. [S ii]/l/l67 17,673 1/Hq' vs. [Nn]/}6584/Ha diagrams. First row: Line ratios for the different LIRGs groups in our sample. 
Second row: Theoretical models used for comparison where those in the mi ddle column have solar metallicity and those in the right 
twice solar metallicity. In all panels, the empirical boundaries proposed in Veilleux & Osterbrockl (Il987h to separate the ionization 
by stars or other mechanisms ass uming a log ([Oiii]/t5007/H/6) - -0.2 are shown. The color code for the models is: i) Green scale: 
Star-formation models from Do pita et al.l (l2006h . Different solid lines represent different R parameter which vary from (darkest 
green) to -6 (lightest green) while dashed lines connect points with th e same age ranging f rom 0.5 Myr (darkest green) to 3.0 Myr 
(lightest green); ii) Orange scale: Dusty AGN models models from iGroves et al.l (l2004l) . Ionization parameter ranges from 0.0 
(bottom) to -4.0 (rightmost edge). Only the cases with a - -2.0 have been plotted. Smaller a parameters typically predict higher 
[S ii]/t/t6717,67 31/Ha and [ N ii]/t6 584/HQ' line ratios for a given ionization parameter, iii) Blue scale: Models for shocks without 
precursor from lAllen et al.l (l2008h for n^ - I cmT^. Solid lines indicate models for magnetic parameters of B/n'^^ = (darkest 
blue), 2, and 4 juG cm-'^^ (lightest blue) while dashed lines joint points with velocities of 100 km s"' (darkest blue), 200 km s"' 
(medium blue) and 300 km s"' (lightest blue). The predictions for models with shock velocities greater than 300 km s"' occupy a 
similar area as those for - 300 km s"' in this diagram. Models for shocks with precursor cover a simila r area in th is diagram 
displaced by ~ -0.2 dex in both line ratios, iv) Black dash-dotted line: Models for turbulent mixing layers of Slavin et al.l (fT99 3) for 
Z - Zq and transverse velocity v, - 25 km s Average temperature for the gas increase with [Nn]/16584/Ha. The three diamonds 
indicates log r=5.0, 5.3 and 5.5. 



Elevated line ratios can in principle be explained by both 
shocks and AGN. This is evident from the models presented in 
Figs. [8] and |9l as wel l as from integrated data from the literature 
dKewlev et al.ll2006l) . Indeed, four of the 32 galaxies are clas- 
sified as Seyfert (the northern member of IRAS F07027-6011, 
IRAS Fl 3229-2934, and the western member of IRAS 14544- 
4255 and IRAS 21453-3511, see footnotes in Table |3]l. The 
northern member of IRAS F07027-6011 should not pose a 
problem, because most of its line ratios are typical of young 
stars. However, as stated in Table |3] some areas with elevated 
line ratios could be caused by the nuclear AGN. The galaxy 
IRAS 13229-3934 contains a central AG N which causes an ion - 
ization cone in the north-south direction teedregal et al.ll2009l) . 
However, our measured line ratios are associated to wide areas 
in the inter-arms zone, not coincident with the direction of the 
ionization cone, and thus difficult to be explained by the cen- 



tral AGN. Finally, for the western member IRAS 14544-4255 
and IRAS 21453-3511, it is not possible to discern with the 
present information whether the ionization in the external areas 
is caused by shocks or the central AGN. However, an indirect 
argument allows us to favor shocks against AGN as the mech- 
anisms responsible of the observed line ratios. As was shown 
for Arp 299, a very nearby LIRG, the extra-nuclear area with an 
excitation caused by the AGN (as seen by the [O ra]/15007/Hj6 
line ratio) is very small and restricted to specific directions de- 
fined b y the ionization con e (see detailed analysis for this sys- 
tem in iGarcfa-Marm et al ] 12006.) . Also, nuclear regions, where 
one could expect the largest influence of an AGN, were removed 
from our analysis (see Sect. 12.41 ). 

It is worth mentioning that the metallicities (solar to twice 
solar) of the models (H ii and shocks) that best fit the range of 
data are underabundant by a factor of ~ 1.5-2.0 with respect 
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Fig. 9. [Oi]/16300/Hq' vs. [Nii]/16584/Hq' diagrams. First row: Line ratios for the different LIRGs groups in our sample. Second 
row: Theoretical models used for comparison where those in the middle column have solar metallicity and those in the right twice 
solar metallicity. The utilized color/symbol code is the same as in Fig. [8] 



to what is expected from the mass-metallicity relation for galax- 
ies dTremonti et alJl2004h . This relation predicts metallicities of 
Z ~ 1.9 - 2.8 2(70 for galaxies with mass es as tho se expected 
in typical LIRGs (~ l O'" - lO'^ Mq, e.g.. lHinz & Rieke 200^ 
IVaisanen et al.ll200 8a b). Similar findings for the nuclea r region 
of the se kinds of galaxies have been reported by .Rupke et alJ 
(I2008h . 

In summary, the ionization of the extended, extra-nuclear re- 
gions in isolated galaxies (class 0), is mostly explained as due 
to young stars like in Hii regions. Systems showing some de- 
gree of interaction (class 1 and 2) present a large and increas- 
ing fraction of regions which are better explained by shocks. 
This is particularly evident when using the the [Oi]/i6300/HQ' 
vs. [Nn]/i6584/HQ' diagnostic diagram. In particular cases (e.g. 
IRAS 21453-3511), an AGN also could explain the observed 
high line ratios. In order to explore the precise relevance of a pu- 
tative AGN, other additional observables than those utilized here 
(e.g. flux in the [Oiii]/I5007 emission line) would be needed. 
Further evidence for the presence and relevance of shocks as 
ionizing sources can be found through the gas velocity disper- 
sions and their correlations with the emission line ratios. This is 
presented in the next section. 



^ W e have employed 12 + log(0/H)o = 8.66, from lAsplund et alj 
ll2004h . 



3.5. The relation between excitation and gas velocity 
dispersion: further evidence for the importance of 
shocks 

The presence and relevance of shocks in (U)LIRGs has already 
been suggested by a positive relation between the ionized gas 
velocity dispersion and its ionization degree as traced by the 

i s n1/l/l67l7,6 73l/HQ' ratio in a sample of about 30 galaxies 
Armus et al.l |l989; Dopita & Sutherland! 119951: l^illeux et al| 
11991 " These studies were based on long-slit observations and 
were therefore dominated by the contribution from the high- 
surface brightness nuclear regions in a large number of cases. 
Also, the slit was positioned along a given orientation and 
therefore the results of these studies do not necessarily rep- 
resent the excitation and kinematics of the ionized gas in the 
extra-nuclear extended (several kpc) regions. Instead, for the 
detailed study of the excitation conditions in the extended re- 
gions it is more appropriate to use the two-dimensional infor- 
mation provided by IFS data once the nuclear regions are re- 
moved, as already discuss ed for a small sample of ULIRGs 
('Monreal-Ibe ro et ani200 6'). We performed a similar study with 
the present sample of LIRGs, limiting the analysis to the data 
with [Nii]/16584/Hq', [S ii]^^67 17,673 l/Ha, and [Oi]A63QQfHa 
higher than -0.7, -0.8, -1.8, respectively (i.e. line ratios in the 
range expected for shocks according to the models discussed in 
the previous section). In particular. Fig. [TO]presents the relation 
between the excitation degree (here represented by our three line 
ratios) and the gas velocity dispersion for the entire sample (up- 
per panels) and according interaction class (three lower panels). 
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Fig. 10. Relation between the velocity dispersion and the [N ii]A65S4/Ua (left), [S ii] /1^671 7,673 l/Ho- (middle), and [O i]/l6300/Ha 
(right). Mean and standard deviation of each 0.2 dex bin in the velocity dispersion are shown with blue crosses and eiTor bars 
respectively. The red dash-3 dotted line in each individual graphic represent the one-degree polynomial fit to these values (see Table 
|4]for specific values). We refer the reader to the On-line only version for the equivalent plots associated with each individual system. 



The relations for each individual pointing have been appended 
in the On-line only version (see Fig. [TT]). 

Several conclusions are already evident from these results. 
Firstly, the ionized gas in the extra-nuclear regions of LlRGs has 
typically velocity dispersions between 32 km s"' (lower limit 
given by our spectral resolution) and 125 km s"', with very few 
regions having velocities above this value, independently of the 
morphology and infrared luminosity. 

Secondly, the velocity dispersion of the ionized gas is larger 
in galaxies with some degree of interaction, i.e. class CI and 
C2. While galaxies identified as isolated have a median veloc- 
ity dispersion of 37 km s"', class CI and C2 have values of 46 
km s"' and 51 km s"', respectively. Moreover, the velocity dis- 
persions of class galaxies (i.e. isolated) tend to be concentrated 



in the low velocity range, while classes 1 and 2 have a relevant 
fraction of regions with high velocities. This means that only 
5% of the data points for class galaxies present velocity dis- 
persions larger than 80 km s while 28% and 16% of those 
for classes 1 and 2 respectively do when considering, for exam- 
ple, the [Oi]/16300/Hq' line ratio. These values are much higher 
than the median velocity dispersions measured in the extranu- 
clear ionized regions o f normal galaxies with velocities in the 
20 to 30 km s" range (lEpinatetalJl2010h . This is a clear indi- 
cation that the ionized interstellar medium in LlRGs in general, 
and even more in interacting LlRGs, is dynamically hotter than 
the quiescent ISM of normal galaxies due to the strong shocks 
produced by the tidal forces and by stellar winds in the powerful 
nuclear starbursts. 
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Thirdly, there is a clear correlation between the excitation 
degree and the velocity dispersion in interacting galaxies (class 
1) and mergers (class 2), while the evidence of correlation in 
isolated galaxies (class 0) is poor and restricted to the relation 
involving the [N ii]/16584/Hq' line ratio. 

To obtain a more quantitative analysis of the degree of 
the correlation, the linear Pearson correlation coefficient, r, is 
used. This coefficient quantifies the degree of coiTelation be- 
tween two given quantities that are assumed to follow a lin- 
ear relation. It varies from -1 to -i-l, where +1(-1) means a 
perfect coiTelation(anti-correlation) and means no-coiTelation. 
The Pearson coefficient for the log([SiiU/16717,6731/HQ-) - 
log(cr) in the 'Armus et al.' (Il989l) sample (see also Fig. 8 in 
[Oopita & Sutherland 1995) has a value of r = 0.4. Hereafter 
we consider that a positive relation exists between the excitation 
conditions and velocity dispersion in the ionized gas only if r 
is higher than 0.4. The computed r coefficient for the different 
groups as well as the polynomial coefficients obtained from a 
least-square one-degree polynomial fit to the mean of the data in 
0.2 dex bins in the velocity dispersion appear in TableH] For the 
individual pointings, the direct fits to the data are included in the 
lower right corner of the corresponding panel in Fig. [17] in the 
On-line only version. 
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Fig. 11. Variation of the r coefficient (up) and percentage of 
galaxies with a confidence level higher than 90% showing corre- 
lation (down) with the interaction types (O=isolated, 1 =interact- 
ing pairs, 2=merger remnants) for the three emission line ratios. 
We used the same line and color code as in Fig.|4]to distinguish 
among the emission line ratios. 



The confidence level of the different relations was esti- 
mated by means of the t-Student test. All relations corre- 
sponding to the different groups present a confidence level 
higher than 99.9%. For individual pointings, the [N ii]/16584/Hq'- 
o" show a larger degree of confidence than those for the 
[Sii]/l/l6717,6731/HQ'-cr and [Oi]^6300/Ha-o- relations, be- 
cause the S/N limits the number of points. All relations in- 
volving the [Nii]/16584/Hq' line ratio except the southern mem- 
ber of IRAS 18093-5744 present a level of confidence higher 
than 95%. For the other line ratios the level of confidence is in 
general lower. For the [S ii]/l/l6717,6731/Ha' line ratio, 11 out 
of 36 pointings (i.e. 30%) did not reach a 90% of confidence 
level. These ai-e IRAS F06295-1735, IRAS Fl 1255-4120, 
IRAS Fl 1506-3851, IRAS FOl 159-4443, ESO 297-G012, 
IC 563, both pointings for IC 564, the northern and central mem- 
ber of IRAS 18093-5744 and IRAS 12596-1529. In the case 
of the relation involving the [Oi]/I6300/Ha' line ratio 22% of 
the pointings (i.e. 6 out of 27) did not reach a 90% confidence 
level. They are the southern member of IRAS F07027-6011, 
ESO 297-GOll, ESO 297-G012, IRAS 08424-3130, the cen- 
tral member of IRAS 18093-5744 and IRAS 12596-1529. 

The variation of the Pearson coefficient (r) with the interac- 
tion class is summarized in the upper part of Fig. [TT] The per- 
centage of galaxies in each group showing a correlation coef- 
ficient r > 0.4 and a level of confidence higher than 90% is 
summarized in the lower part of Fig.[TT] The quantitative analy- 
sis emphasizes the main differences according to the interaction 
class already mentioned above. The degree of correlation be- 
tween the excitation properties and the velocity dispersion of the 
ionized gas when the [S ii]/l/167 17,673 1/Ha and [Oi]A63QQfHa 
line ratios are considered increases with the degree of interac- 
tion. In particular, the [Oi]/I6300/Ha' ratio is the best tracer of 
shocks as already shown in the previous section. When using this 
ratio, the Pearson correlation coefficient for the class galaxies 
is r = 0.30, indicating no correlation. This result supports the 
idea that shocks are not playing an important role in the ioniza- 
tion of the external areas of this group of LIRGs and agrees well 
with our previous findings that young stars are the main cause for 
the ionization. On the other hand, the mean r values for classes 1 
and 2 (0.48 and 0.61, respectively) are well above the 0.4 crite- 
rion, which we considered necessary for a good positive relation. 
This quantitative result clearly supports a direct cause-effect re- 
lation between the dynamical status of the gas, i.e. turbulence 
and shocks traced by the velocity dispersion, and its excitation 
conditions. 

In summary, there are two clear differences between LIRGs 
classified as class (i.e. isolated), and those classified as class 1 
(mostly pairs), and class 2 (mergers). The ionized gas in classes 
1 and 2 is characterized by covering a wide range in velocity 
dispersion with an extension towards higher values than class 
galaxies. Moreover, the dynamical status of the gas, turbulence 
and shocks, plays an important role in the excitation of the gas 
mainly in LIRGs classified as interacting pairs or evolved merg- 
ers. 

3.6. Tidal forces as tine origin of sliocks in fhe extra-nuclear 
regions? 

It was shown above that the importance of shocks in large ex- 
tended regions increases with the interaction class. What is the 
origin of these shocks? Are the shocks caused by the tidal forces 
due to the interaction process itself? Or are the shocks produced 
in stellar superwinds associated with the intense starbursts gen- 
erated in the nuclear regions? 
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Fig. 12. Same as Fig-fTOlbut for class 1 and 2 systems grouped in four luminosity bins. 



To answer these questions, we investigated how the rela- 
tion between the excitation conditions and gas velocity disper- 
sion varies with the total infrared luminosity of the system. In 
the sample considered here, only four objects were classified as 
Seyfert - and thus h arboring an AGN- according to their nuclear 
specti-a ( Kewlev et al. 2001b; Corbett et al. 2003) (see Table|3]l. 
Thus the total infrared luminosity is considered in general a di- 
rect tracer of the intensity of the star formation, because the in- 
frared luminosity i n these galaxies s cales linearly with the star- 
formation rate (e.g. Kennicutt 1998). Figure [12] presents similar 
relations to those in Fig. [10] for those classes with some degree 
of interaction (class 1 and 2) but this time binned in four lumi- 
nosity ranges covering the entire luminosity range of LIRGs (i.e. 
from 10" to 10'2 Lq). 

Contrary to what happens with the interaction class there 
is no evidence for a dependence of the correlation coefficient 
with the luminosity bin, independently of the utilized line ratio. 



Assuming star-formation dominates the energy output in these 
galaxies (see above), the range in infrared luminosity covered 
by the sample repr esents a change of a factor ten in the star- 
formation rate (see iKennicuttI [19841 for the specific relation). 
Therefore the radiative and mechanical energy released in the 
surrounding ISM due to supernovae explosions and stellar winds 
produced in young massive star s would increase lin early with the 
rate of star formation (see e.g. 'Colin a et al.lll99 1^. These linear 
relations in the mass, momentum and energy deposition have ai- 
re adybeenjneasuredjn the cool, neutral gas traced by the Na i 
line (iRupke et al.ll2005bl) . However, these studies do show some 
evidence for a flattening of these relationships for star-formation 
rates above 10 Mq yr"' (i.e. LIRG and ULIRG range). It is un- 
clear whether these results would apply to the ionized gas traced 
by the Ha line. The momentum and energy release in the cool, 
neutral gas is usually a small fraction of the energy in the warm 
ionized and hot X-ray emitting gas. Moreover, the nonlinear re- 
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lation betw een the outflo w velocity of the cool gas and the SFR 
(ocSFR°^^ 'Martinl l2005h could indicate a saturation in the me- 
chanical energy liberated into the more dense neutral gas, but not 
necessarily in the other phases of the ISM. In addition, the veloc- 
ity gradients measured in several outflows detected in ULIRGs 
are inconsistent with the expected gaseous radial flows pr oduced 
by nuclear (size of 200-300 pc) starbursts ('Ma rtinll2006l) . These 
gradients could indeed still be consistent with more extended 
starbursts on scales of kpc, or s hocks generat ed by tidal forces 
during the interaction process dMartinl 120061) . Our Ha IFS of 
(U)LIRGs indicates that the Ha emission is more concentrated 
than that of the stellar continuum. In particular, the fraction of 
Ha emission within the central 2 kpc is higher than that of the 
continuum for about 80% of the cases. However, ~6Q% of the 
objects have more than half of their Ha emission outside the 
central 2 kpc (see Rodriguez-Zaurm et al. in prep, for details). 

Therefore, if the detected evidence of shocks were due to 
star formation, it would be reasonable to expect a more turbu- 
lent gas with stronger outflows and shocks. This could be traced 
by higher excitation conditions and velocity dispersions in the 
ionized gas as well as higher correlation degrees. No evidence 
for any of this is observed in the present data that sample large 
extra-nuclear regions of several hundreds of pc to several kpc in 
size, outside the circumnuclear regions. 

The comparison between Figs. [TOl and fT2] indicates that the 
presence and relevance of shocks are more strongly coiTelated 
with the interaction/merging class of a system than with its 
star formation activity. Indeed, our sample of class galaxies 
shows on average a level of star formation activity (i.e. infrared 
luminosity) similar to or slightly lower than that of classes 1 
and 2 (log(L/Lo) = 11.44 for class against 11.51 and 11.60 
for classes 1 and 2, respectively). Note that this statement is 
valid for our particular sample, which was selected in order to 
cover all the interaction types and luminosity ranges in a more 
or less uniform way. However, it does not apply to complete 
samples of LIRGs, because they show a much higher percent- 
age of interaction /merging systems at higher luminosities (e.g. 
ISanders & Ishidal i2004). These results point to tidal forces asso- 
ciated with the interaction/merging process as the origin for the 
shock ionization in the extended, extra-nuclear regions. A sim- 
ilar result was found by MAC06, who concluded that the more 
likely explanation for shocks in five out of the six ULIRGs stud- 
ied there were tidally induced large scale gas flows caused by 
the merging process. 

D etailed studies of the ne arest ULIRG Arp 220 dColina et al.l 

120031; iMcDoweU et al.ll2003l) have also suggested that large ex- 
tended regions in this system arise purely from merger dynam- 
ics and collisional shock heating of the gas. At the same time, 
footprints associated wit h starburst s uperwinds ha ve also been 
detected (iHeckman etaDflgQO- . AiTibas et al.ll200lh . In Arp 220, 
the ionized gas plumes that could be associated with superwinds 
generated in the nuclear starburst form an elongated structure up 
to a distance of about 2 kpc from the nucleus. However, for the 
present sample of LIRGs, which show a significant lower star- 
forming activity based on their L/^ but sirti ilar dynamical mass 
(iHinz & Riekell2006HVaisanen et alJl2008bllal) and hence similar 
escape velocity, a smaller area of influence of the stellar super- 
winds should be expected. 

Thus our results indicate that the tidal forces during the inter- 
action process are the mechanism producing the ionizing shocks 
in the extended extra-nuclear regions in LIRGs. This is still com- 
patible with the existence of SGWs produced in nuclear star- 
bursts at distances closer to the nucleus (i.e. radius of 1 to 2 



kpc), or AGN ionizing cones alon g particular orientations (e.g. 
Arp 299. lGarcia-Marin et al.ll2006l) . 

3.7. Interacting LIRGs and ULIRGs. Towards a common 
logCP \]A6300/Ha) - log(cr) relation? 

In previous sections we showed that our IFS data indicate that 
shocks produced by the tidal forces in interacting and merging 
LIRGs play a relevant role in the excitation of the extended ion- 
ized regions, without any clear relation with the intensity of the 
star formation. Because ULIRGs are the extreme cases of inter- 
actions and mergers, one may wonder how the LIRGs situation 
compares with that for the ULIRGs. For that purpose, we com- 
pared our results with those presented in MAC06. The relatively 
small sample of MAC06 (i.e. six ULIRGs systems, nine galax- 
ies) was made out of three class 1 systems (IRAS 08572+3915, 
IRAS 12112+0305, IRAS 14348-1447) and three class 
2 systems (IRAS 15206+3342, IRAS 15250+3609, 
IRAS 17208-0014). We found that in afl the systems but 
IRAS 17208-0014 the extra-nuclear extended regions were 
well explained by ionization due to fast shocks with velocities of 
150 - 500 km s Because none of the systems in the MAC06's 
ULIRGs sample is classified as isolated (class 0), we present in 
Fig. [T3] (second row) the line ratio vs. velocity dispersion only 
for LIRGS of classes 1 and 2. 

Considering the best shock tracer, the [Oi]/I6300/Ha line 
ratio, the similarity of the log([Oi]/l6300/Ha) - log(cr) correla- 
tion found independently for LIRGs and ULIRGs is remarkable. 
Because the average infrared luminosity (\og{Lir I Lq)) of the 
ULIRG sample is 12.24, this extends the previous result found 
for class 1 and class 2 LIRGs, where a similar linear relation 
exists for the entire LIRG luminosity range, independent of the 
luminosity beam selected. Although the number of ULIRGs in 
the MAC06 sample is small, the combined results suggest a com- 
mon relation between the excitation and dynamical properties of 
the ionized gas in interacting and merging (U)LIRGs over the en- 
tire infrared luminosity range of log(L,>/LQ) - 11.1 - 12.3 cov- 
ered by these two samples. This relation is best traced by the di- 
rect proportionality (slope of ~ 1 .0) between log([0 i]/16300/Ha) 
and log(cr) (see Tables |4] and |6] for the specific slopes). A larger 
sample of about 20 ULIRGs with available optical integral field 
spectroscopy is currently studied (Garcia-Marm et al. in prepa- 
ration) to further investigate the reality of the relation for both 
LIRGs and ULIRGs and its physical interpretation. 



4. Conclusions 

The two-dimensional ionization structure of the extended (few 
to several kpc) ionized gas in a representative sample of 32 low- 
z LIRGs (i.e. log{Li,-/LQ) = 11 .00 - 12.00 luminosity range) was 
investigated with the VIMOS integral field spectrograph. The 
sample covers isolated galaxies, as well as interacting galaxies 
and systems in an advanced stage of the merger. This paper in- 
vestigates the nature and origin of the main ionization mech- 
anisms operating in the extra-nuclear regions of these systems 
based on several thousands of independent measurements of the 
emission line ratios (up to twenty 24000 for [Nii]/i6584/Ha) 
and velocity dispersions. 

The present study is part of a larger project devoted to the 
study of the two-dimensional structure for the stars and ionized 
gas as well as its kinematics and ionization conditions in repre- 
sentative samples of LIRGs and ULIRGs using optical IFS. The 
main results of this study can be summarized as follows: 
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Fig. 13. Upper row: Same as Fig. [TO] but for all systems in MAC06 but IRAS 17208-0014. This is similar to the first row of Fig. 
4 of MAC06 but we plotted the same ranges in the axes as in Fig.[TO]for a better comparison with the results in the present paper. 
Lower row: Same as Fig.[TO]for the combined class 1 + class 2 group. 



1. The distribution of the [Nii]/16584/Ha line ratio does 
not show any significant variations with the interaction class, 
with most regions presenting a ratio typical of Hii regions. The 
[S ii]i^67 17,673 1/Hq' and [Oi]A63Q0fHa line ratios do however 
show a change in their distribution with an extension towards 
higher excitation (i.e. LINER-like excitation) for galaxies clas- 
sified as interacting pairs and advanced mergers. This change is 
more pronounced for the [O i]A63QQfHa ratio. 

2. There is an anti-coiTelation between the ionization degree 
and the Ho' surface brightness independently of the interaction 
type and similar to what occurs in our Galaxy or in the so-called 
DIG in other spiral galaxies. Most of the observed line ratios are 
similar to those found in Hii regions in our Galaxy, but there 
is a relatively large percentage of fine ratios similar to those for 
DIGS. 

3. The nature of the ionization sources was investigated com- 
paring the measured [S ii]/1/16717,6731/Hq' vs. [Nii]/16584/Hq' 
and [Oi]/16300/Ha vs. [Nn]/i6584/HQ; line ratios with the pre- 
dictions of ionization due to stars (Hii regions), TML, shocks 
(DIGs) and power-law (AGN) spectra. Turbulent Mixing Layers 
do not seem to play a major role in the ionization of the extra- 
nuclear regions. Line ratios in LIRGs classified as isolated can 
mostly be explained as caused by ionization due to young stars. 
On the other hand, the ionization in a large fraction of the re- 
gions in systems with some degree of interaction cannot be due 
to stars but is better explained by high velocity shocks. This 
is particularly evident when using the best shock tracer, i.e. 
the [Oi].16300/Hq' vs. [Nn]/}6584/Ha diagram. Independently 
of the ionization mechanisms, only models with metallicity be- 
tween solar and twice solar are able to explain the observed line 
ratios. 

4. Local velocity dispersions increase with the interaction 
degree, with medians of 37, 46, and 51 km for class 0, 1, 
and 2 respectively, and are highe r than those t ypically found in 
normal spirals (~ 20 - 30 km s '. lEpinat et al.l 2010). This indi- 
cates that the ionized ISM in LIRGs is dynamically hotter than 



the quiescent ISM of normal galaxies due to strong shocks pro- 
duced by tidal forces and stellar winds associated to the nuclear 
starburst. 

5. There is a positive relation between the degree of exci- 
tation (as traced by the emission line ratios) and the velocity 
dispersion of the ionized gas in LIRGs classified as interact- 
ing systems, and mergers, while this relation is not observed 
in isolated systems. This relation is better seen when using the 
[Oi]^6300/Hq' and [S ii]^^6717,6731/Ha line ratios, and sup- 
ports the scenario where the relevance of shocks as ionizing 
sources in the extranuclear extended regions of LIRGs increases 
when there is some degree of interaction. 

6. The relation between the degree of excitation and the ve- 
locity dispersion of the ionized gas in interacting and merging 
LIRGs does not clearly improve with the infrared luminosity (i.e. 
star formation rate) of the systems. Thus the interaction process 
itself rather than superwinds caused by the star formation seems 
to be the main origin of the shocks in the extended extra-nuclear 
regions, assuming that the release of energy into the ISM is pro- 
portional to the SFR. This result is still compatible with stellar 
superwinds in the internal regions of these systems, and/or along 
certain preferential directions associated with AGN-related out- 
flows. 

7. A comparison between the sub-sample of interact- 
ing/merging LIRGs and a small sample of ULIRGs suggests the 
existence of a common positive log([Oi]/l6300/HQ') - log(cr) re- 
lation. If confirmed, these results will provide further evidence 
for the tidal origin of shocks in these galaxies over the entire 
LIRG and ULIRG luminosity range. A study with a larger sam- 
ple of ULIRGs is under way to confirm the result. 
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Table 1. VIMOS LIRG sample: For each pointing, the total number of data points, mean value and standard deviation for the 
different Une ratios and velocity dispersion are given. 



Galaxy 




D 


Scale 


[Nii]/16584/Hq' 


[Sii]ii6717,6731/HQ' 


[0 i]^6300/Hq' 




(kms-') 


(IRAS number) 


(Le) 


(Mpc) 


(pc/") 


n 


X ± cr(x) 


n 


X ± cr(x) 


n 


X ± cr(x) 


n 


~ + cr(x) 


Class 


F06295-1735 


11.27 


92.7 


431 


998 


0.38±0.07 


486 


0.31 ±0.06 


194 


0.05+0.04 


1118 


33+11 


F06592-6313 


11.91 


100.0 


464 


199 


0.82±0.32 










266 


54+33 


F07027-6011 N 


11.64 


137.4 


626 


342 


0.40+0.09 


94 


0.30+0.06 


92 


0.07+0.04 


387 


32+11 


F07027-6011 S 




137.4 


626 


163 


0.61±0.17 


74 


0.24+0.05 


49 


0.03+0.02 


177 


72+17 


F07160-6215 


11.16 


46.7 


221 


755 


0.92±0.61 


189 


0.33+0.16 






888 


69+37 


F10015-0614 


11.77 


73.1 


343 


1298 


0.41±0.10 


850 


0.31+0.08 


137 


0.03+0.01 


1312 


51+13 


F10409-4556 


11.26 


91.4 


425 


928 


0.49±0.25 










1083 


39+16 


F10567-4310 


11.07 


74.6 


350 


1154 


0.44±0.15 


270 


0.25+0.06 






1153 


29+11 


Fl 1255-4120 


11.04 


70.9 


333 
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0.28+0.06 
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33+20 


F11506-3851 


11.30 


46.6 
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0.67±0.23 
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0.33+0.13 
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41 + 12 
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61 


0.02+0.01 
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Class 2 
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11.69 


69.1 
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113.1 
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11.77 


149.9 
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0.13+0.09 


303 
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61+19 


12116-5615<3' 
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0.78±0.28 


80 
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70.0 


329 
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43+19 



*'* Logarithm of the infrared l uminosity, L,> = L(8 - lOOOyum), in units of solar bolometric luminosity, calculated using the fluxed in the four 
IRAS bands as given in Sander s et al.l ( l2003h . when available. Otherwise from the IRAS fluxes given in IRAS Point Source and Faint Source 
catalogs Mosbir & et al. ( 1990). When a system of more than one galaxy is observed using several VIMOS pointings, the infrared luminosity of 
the entire system is indicated in the data corresponding to the first galaxy of the system. 
Velocity dispersion derived from the Hor line. 

Classification uncertain. Not used in the analysis of systems grouped by interaction class. 
Members of IRAS 09437+0317. 



Table 2. Median, standard deviation, and number of data points for the distributions of the different groups and line ratios. 



Group 


log([Nn]^6584/H(!r) 


log([S 


ii]/l/l67 17,673 1AI«) 


log([Oi]^6300/Ha) 




Median 


Std Dev 


N 


Median 


Std Dev 


N 


Median 


Std Dev N 


All 


-0.37 


0.18 


24017 


-0.49 


0.16 


12 564 


-1.46 


0.34 5 733 


Class 


-0.37 


0.18 


10 605 


-0.55 


0.13 


4932 


-1.56 


0.28 1 953 


Class 1 


-0.38 


0.17 


6680 


-0.49 


0.18 


3 140 


-1.43 


0.33 1 349 


Class 2 


-0.36 


0.16 


3 992 


-0.48 


0.17 


2 628 


-1.35 


0.35 2 023 
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Table 3. Ionization in the extra-nuclear ionized regions of LIRGs. 
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Galaxy 


Mechanism"'* 




Comment 


(IRAS number) 








Cla,ss 


rUOZyj— I /JJ 


Young stars 


0.22 




rUDJyZ— UJ 1 J 




0.42 


iNo [u ij/tujUu nor [o iij/t/(o / 1 / ,D/ 31 data, r based on [JN nj/iojfi4/riQ'. 


rU/UZ/— DUil In 


Young stars and snocks ( 1 00-200 km s ) 


44 


Part of the shocks may be associated to the nucleus according [0 i]/l6300/Ho' map. 


F07027-60I1 S 


Young stars 


0.21 


r based on [S u]AA67 1 7,673 l/Ha. r based on [0 i]i6300/HQ' not reliable. 


F07 160-62 15 


Shocks (100-300 km s"') 


0.31 


r based on [S n]i.l67 1 7,673 l/Ha. No [0 1]^6300 data. 


F100I5-0614 


Young stars 


0.15 




F10409-4556 




0.37 


r based on [N ii]^6584/Ho'. No 10 i]i6300 nor [S ii] A Abl 11, 61 3 1 data. 


F10567-4310 


Young stars 


0.23 


/■ based on [S n]^i67 1 7,673 l/Ha. No [0 1]^6300 data. 


Fl 1254-4120 


Young stars 


0.71 


r based on IN n]i6584/HQ'. r based on [S n]/l^67 17,6731/00- not reliable. 








No 10 r]yl6300 data. Possible shocks associated to the bar according to the 








[S u]AA61 1 7,673 l/Ha map. 


Fl 1506-3851 


Young stars 


0.19 


Some small areas better explained with low velocity shocks. 


F121 15-4656 


Young stars 


-0.24 




F13229-2934"" 


Shocks (100-300 km s"') 


0.47 


Some areas explained by stars, mostly associated with the aims. Part of the high 








excitation may be associated to the nucleus according to the [0 i]/i63()0/Hti' map. 


F22 132-3705 


Young stars 


0.10 


Small areas better explained with low velocity (<100 km s"') shocks. 


Cla,ss 1 


FO 1159-4443 


Low velocity (< 100 km s ') shocks 


0.23 


Some regions explained by stars. 


ESQ 297-GOll 


Young stars 


-0.12 


Some small areas explained by shocks, r based on [S n]/l/1671 7,673 l/Ho-. r based 








on [Oi]/1630()/Ho' not reliable. 


ESQ 297-G012 


Young stars and shocks 


-0.56 


r based on [N n],16584/H£i'. /■ based on [0 i],16300/H£i' or [S n],1.167 17,673 l/Ha 








not reliable. 


F06076-2139 




0.75 


r based on [N n]i6584/HQ' data. No [0 i]/i63()0 nor [S n].li67 17,673 1 data. 


F06259-4708 PI 


Young stars 


0.74 


Low velocity ( 1 00- 1 50 km s~ ' ) shocks also present. 


F06259-4708 P2 


Shocks (100-150 km s"^) and young stars 


0.31 




08424-3130 


Stars and shocks (~ 100 km s"') 


0.61 


r based on [S n]^^67 17,673 l/Ha. r based on 10 i].^6300/Hq' not reliable. Northern 








galaxy mainly shocks; southern one mainly stars. 


F08520-6850 


Stars and shocks (~ 100 Ion s"^) 


0.40 




IC563 


Young stars 


0.22 


Low velocity (100-150 km s~' ) shocks could be also present. 


IC564 N 


Young stars 


0.44 


Low velocity {100-150 km s~') shocks also present, r based on lNn]/16584/H(y. 








No lOi]/16300 data, r based on [S n]^.1671 7,673 l/Ho- not reliable. 


IC564 S 


Young stars 


0.34 


r based on [N n],16584/HQ-. No [0 1],16300 data, r based on [S n],1.167 17,673 l/Ha 








not reliable. 


12042-3140 


Young stars and shocks (100-150 km s"^) 


0.81 


Northern galaxy mainly stars. Southern one shocks. 


12596-1529 


Young stars and low velocity shocks 


0.47 


r based on [S n]^.167 17,673 l/Ho- data. No [Oi]^6300 data, r based on 








[S n],1.167 17,673 l/Ha not reliable. 


F14544-4255 E 


Young stars 


0.67 


r based on [S n]/U67 1 7,673 l/Ha. No [0 1]^6300 data. 


F14544-4255 W*"" 


High velocity (100-300 km s"') shocks 


0.67 


Shocks may be associated to the nucleus according to the [0 i]/i6300/HQ' map. 


18093-5744 N 


Young stars 


0.09 


Low velocity shocks also possible. 


18093-5744 C 


Stars older than 3 Myr and/or with Z < Z© 


-0.65 


rbased on [S n]^i671 7,673 l/Ho-. rbased on 10 i]^6300/H(i' not reliable. 


18093-5744 S 


Young stars 


0.27 




Class 2 


F043 15-0840 


Shocks (100-200 km s ' ) and young stars 


0.63 


Young stars in the inner region and arms/tidal tails. Shocks at 2-3 kpc from the 








nucleus. 


08355-4944 


Young stars 


0.14 


Low velocity shocks also present at ~ 2 kpc from the nucleus. 


F10038-3338 


Shocks (100-200 km s"') 


0.50 


Shocks may be associated to the nucleus. 


10257-4338 


Shocks (100-150 km s"' ) and stars 


0.63 




12116-5615 


Young stars and low velocity shocks 


0.57 




13001-2339 


Shocks (100-150 km s"') 


0.48 


All data with [0 i]/16300/Hq'> -1.2. Part of the shocks may be associated to the 








nucleus 


17138-1017 


Young stars and low velocity shocks 


-0.23 


No tOi]i6300data; r according to [S n]^,16717,6731/Hff data 


21453-3511"" 


High velocity (100-300 km s"') shocks and/or AGN. 


0.76 


Young stars (spiral arms). Shocks may be associated to the nucleus according 



[O i]/i6300/Hg map. 



At least otherwise indicated, based on the comparison of the data with several ionization models in the [O i]/16300/Hq-- [N n],16584/Ha diagnostic diagram. For the present analysis solar 
metallicity has been considered as baseline. If higher metallicity (e.g. Z— 2 Z©) is considered, a higher presence of shocks with respect to stars would have been obtained (see Figurel9l. 

Pearson coefficient inferred from the [O i]-I6300/Hq'- it relation (see text). For objects where [O i]/16300 was not detected or the r coefficient did not reach a confidence level of 90%, r 
was obtained from the [S ii]/(,I67 1 7.673 1/Hff or [N n]/16584/Hff, as indicated. 
"1 Classified as Seyfeit ( Kewley et al. 2001a). 

Classified as Seyfeit tCorbett et aU2003i) . 



Table 4. Linear Pearson correlation coefficients r and 1 -degree polynomial fits for the different morphological groups considered 
here. 



Group 


[Nii]/16584/Hq' 


[Sn]^67 17,673 I/Hq- 


[O i]A6300/lla 




r 


A-l-Bx 


r 


A-i-Bx 


I 


A-t-Bx 


AH 


0.43 


-1.08 -H 0.44 X 


0.33 


-1.44 -1- 0.55 X 


0.54 


-3.12-1- 1.02 X 


Class 


0.54 


-1.17 + 0.54 X 


0.23 


-1.00-1- 0.28 X 


0.30 


-2.95 -1- 0.88 X 


Class 1 


0.43 


-1.23 0.50 X 


0.33 


-1.55 -1-0.60 x 


0.48 


-3.14-1- 1.01 X 


Class 2 


0.31 


-0.79 -1- 0.27 X 


0.60 


-1.18 -H 0.69 X 


0.61 


-3.00 + 0.97 X 
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Table 5. Linear Pearson correlation coefficients r and linear fits for for the class 1 and 2 LIRGs distributed in luminosity groups 
spanning the entire luminosity range. 

Luminosity Range [N ii]^6584/Hf>' [S ii]i67 17,673 l/Ho- [Oi]i6300/Ha 

(log(L/Lo)) r A+Bx r A+Bx r A+Bx 

11.00-11.25 049 -0.97 + 0.39 x 063 -1.39 -f 0.54 x ^005 -1.12 -O09x 
11.25-11.50 0.72 -1.24 4- 0.51 x 063 -1.19-1- 037 x 068 -3.09 1.01 x 
11.50-11.75 0.28 -0.77 -I- 0.22 x 0.24 -081 -I- 0.20 x 035 -2.47 -I- 0.61 x 
11.75-12.00 0.59 -1.58 H- 0.64 X 0.43 -1.21 -h 0.41 x 0.46 -3.01 h- 0.95 x 

Table 6. Linear Pearson correlation coefficients r and linear fits for all CI and C2 LIRGs and ULIRGs 



Group [N II] i65 84/^1- [S ii]^67 17,673 1/Ha [Oi]i6300/Ha 

r A+Bx r A-l-Bx r A-l-Bx 

ULIRGs MAC06 061 -1.99 -i- 082 x 051 -1.36 0.50 x OTTl -3.15 -i- 1.11 x 

LIRGs C1-I-C2 0.38 -0.96 -I- 0.35 x 0.41 -1.09 -i- 0.36 x 0.45 -2.74 -i- 0.79 x 
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FULL VERSION at |http://vrww.damir.iem.csic.e s/extragalactic/publications/p ublications.html| 

Fig. 14. Panel showing the images utilized for the morphological classification. The left column contains details about the size and 
scale for the displayed images, as well as relevant morphological features and final classification. The central and right columns 
display the Digital Sky Survey and HST images, respectively. Orientation is north up, east to the left. 
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Fig.ll4l continued. 
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Fig.ll4l continued. 
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Fig.ll4l continued. 
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Fig.ll4l continued. 
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Fig.ll4l continued. 
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Fig.ll4l continued. 
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Fig. 15, [S n]A/l6717,6731/Ha vs. [N ii]/l6584/Ha diagrams for the individual pointings. 
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Fig.1151 continued. 
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Fig. 16. [O i]/16300/Hq; vs. [N ii]/16584/Hq' diagrams for the individual pointings. 
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Fig.ll6l continued. 
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Fig. 17. Line ratios vs. velocity dispersions relation for the individual pointings. 
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Fig.1171 continued. 
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Fig.1171 continued. 

FULL VERSION at h ttp://www.dami r.iem.csic.es/extragala ctic/publications/publications.htmll 

Fig.1171 continued. 
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Fig.1171 continued. 
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Fig.1171 continued. 
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Fig-HZl continued. 



